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Abstract. The cosmological evolution in the radiation dominated regimen is usually computed 
by assuming an ideal relativistic thermal bath. In this note, we discuss the deviation from the non- 
interaction assumption. In either the standard model (SM) and the minimal supersymmetric standard 
model (MSSM), the main contribution comes from the strong interaction. An understanding of these 
effects are important for precision measurements and for the evolution of scalar modes, where the 
commented corrections constitute the main source of the dynamics. 
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THE RELATIVISTIC THERMAL BATH 

We can describe the thermodynamical properties of any fluid by its free energy den- 
sity (J^). Equivalently, we can use gf, the effective free energy number of relativistic 
degrees of freedom, that is defined by normalizing to the free energy density of one non- 
interacting massless bosonic particle [HI]. The partition function and gf can be computed 
diagrammatically [2]. The one loop vacuum diagrams provide the free theory result. In 
this case, the coefficient gf h simply the number of bosonic {Nb) and fermionic {Nf) 
relativistic degrees of freedom: gfj^^ = Nt, + 7Nf/S, where the 7/8 coefficient multi- 
plying the fermionic contribution is due to the different value of the partition function 
obtained from the Fermi-Dirac thermal distribution rather than the Bose-Einstein one. 
Higher loops account for the interactions. In [1], it has been shown that the strong in- 
teraction is dominant even at very high temperatures (see Fig. [1]) except for the case of 
a heavy Higgs boson. Unfortunately, the situation is not simple since the finite temper- 
ature perturbation approach to QCD has a very slow convergence [3]. In any case, the 
leading order correction coming from the strong interaction has the right order and it can 
be used to estimate the interaction effects in the thermal bath (read yj for a particular 
discussion inside the SM and the MSSM) 

In the right panel of Fig. [T] we show different contributions to the departure of the 
equation of state from the w = 1/3 value corresponding to the non-interacting case. 
The effect decreases as the temperature (T) increases, since 053 decreases with T (and, 
consequently the absolute value of the two loop contribution to gf decreases). At the 
lowest temperatures shown, the departure is stronger in the SM than in the MSSM 
case. This is because as runs faster in the SM. However, since as decreases less in 
the MSSM, the departure from the non-interacting value becomes greater for the MSSM 
as the temperature increases (see Fig.[T]and yj]). 
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FIGURE 1. Left panel: The absolute value of the leading corrections (two-loops) to the effective free 
energy number of relativistic degrees of freedom from the different gauge interactions, top Yukawa and 
Higgs self-coupling in the SM. The strong contribution (0:3), is dominant and is also plotted for the case 
of the MSSM assuming that all supersymmetric partners have masses fixed at otsusy = 500 GeV. All 

(2) 

contributions to are negative. Right panel: The absolute value of the leading correction to 1/3 — CO 
from the same interactions shown in the left panel. The strong contribution is dominant not only due to 
the fact that the correction to gj: is bigger, but also because the running of is steeper. On the other 
hand, the contribution from the hypercharge interaction (and the Higgs coupling for high temperatures) is 
negative and opposite to the others since ai (X) increases with T. Read UJ for more details. 



Precision Cosmology 

Interaction effects in the relativistic thermal bath are important in the precision era 
of astrophysical observations that we have already entered. We can illustrate this ef- 
fect with the thermal relic abundance. Although there are other possibilities [4], Dark 
Matter (DM) is usually assumed to be in the form of stable Weakly-interacting massive 
particles (WIMPs) that naturally freeze-out with the right thermal abundance. WIMPs 
emerge in different well-motivated particle physics scenarios as in R-parity conserv- 
ing supersymmetry (SUSY) models 11,0], universal extra dimensions (UED) IvLS], or 



brane-worlds 



The present uncertainty on the total DM density using data from the Wilkinson Mi 



crowave Anisotropy Probe (WMAP) is around 6% (five years data 11211 ). This precision 
can be improved to 3% when measurements of Baryon Acoustic Oscillations and Type 
la supemovae are also taking into account [13]. On the other hand, the prospects for the 
precision that can be obtained by analyzing data from the Planck Surveyor, reduce this 
uncertainty even below 1%. However, in order to associate the observed precision with 
the thermal relic density of a particular WIMP model, it is neccesary to understand the 
thermal bath up to the same level. Unfortunately, the uncertainties are bigger than 1% 
percent even at high temperatures. The situation is even worse because typical WIMPs 
freeze out at a temperature of 10 GeV approximately, i.e. much closer to the QCD phase 
transition. In this case, it is clear that the perturbative approach does not converge and 
non-perturbative analises (such as lattice studies) are necessary. Alternatively, if we are 



able to understand the nature of DM from other experiments as the new generation of 
colliders [14], we may have the oportunity to improve our knowledge about QCD by 
making precise astrophysical observations. 



Scalar Dynamics 

On the other hand, there is a case in which the deviation from the ideal gas is 
not only important but also fundamental, since it constitutes the leading order term 
of the dynamics. It happens for the case of scalar fields that couple to the trace of 
the energy momentum tensor. These scalars are ubiquitous in many theories, as for 
example, Jordan-Fierz-Brans-Dicke (JFBD) scalars in Scalar-tensor theories of gravity 
such as f(R) theories [|l5|l . moduli fields from string theories such as the dilaton [16], or 
graviscalars from extra dimensional models such as the radion (v^. 

The dynamics of these theories is usually described in what is called the Jordan Frame 
or Einstein Frame. Inside the Jordan Frame, the standard matter is coupled minimally 
to the Jordan metric. It implies that the particle masses are constant, whereas the Planck 
scale depends on the particular value of the scalar field. Alternatively, it is possible to 
use the Einstein Frame, where the action for the metric is the standard, i.e. it is described 
by the Einstein-Hilbert action, but the matter content is coupled explicitly to the scalar 
mode. In this frame, the Planck scale is constant but the particle masses depend on the 
value of the scalar field. 

Both descriptions are classically equivalent and are related by a simple conformal 
transformation. The cosmological evolution of the scalar field is determined by the 
standard equation of motion of a scalar field with an extra source term proportional 
to the trace of the energy momentum tensor. If the relativistic thermal bath is ideal, this 
term is zero. It means that the scalar field is fundamentally frozen at its initial value. 

The situation is completely different if this source term is non-zero, even if it is 
small. The conformal factor may evolve orders of magnitude depending on the particular 
coupling. For example, in [1], it has been studied a model with quadratic coupling and 
without potential. In this case, the coupling works as an attractor to the general relativity 
limit. This important effect implies that the early cosmology can be very different from 
standard. For example, depending of initial conditions, there is a maximum in the Jordan 
temperature that can solve or alleviate the unwanted relic problem since it could avoid 
the production of dangerous solitons or unstable gravitinos [IJ. 

In conclusion, we have discussed briefly the deviations from the ideal relativistic ther- 
mal bath, with interesting consequences for precision measurements and the dynamics 
of scalar modes. 
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